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On the reachability and observabihty 
of path and cycle graphs 

Gianfranco Parlangeli Giuseppe Notarstefano 

Abstract 

In this paper we investigate the reachabihty and observabihty properties of a network system, 
running a Laplacian based average consensus algorithm, when the communication graph is a path or 
a cycle. More in detail, we provide necessary and sufficient conditions, based on simple algebraic 
rules from number theory, to characterize all and only the nodes from which the network system is 
reachable (respectively observable). Interesting immediate corollaries of our results are: (i) a path graph 
is reachable (observable) from any single node if and only if the number of nodes of the graph is a 
power of two, n = 2% i e N, and (ii) a cycle is reachable (observable) from any pair of nodes if and 
only if n is a prime number. For any set of control (observation) nodes, we provide a closed form 
expression for the (unreachable) unobservable eigenvalues and for the eigenvectors of the (unreachable) 
unobservable subsystem. 

I. Introduction 

Distributed computation in network control systems has received great attention in the last 
years. One of the most studied problems is the consensus problem. Given a network of processors, 
the task of reaching consensus consists of computing a common desired value by performing 
local computation and exchanging local information. A variety of distributed algorithms for 
diverse system dynamics and consensus objectives has been proposed in the literature. 

An early short version of this work appeared as yj: differences between this early short version and the current article include 
the reachability analysis, a much improved comprehensive and thorough treatment, revised complete proofs for all statements. 

The research leading to these results has received funding from the European Community's Seventh Framework Programme 
(FP7/2007-2013) under grant agreement no. 224428 (CHAT) and n. 231378 (C03AUV) and from the national project "Sviluppo 
di nuovi metodi e algoritmi per I'identificazione, la stima bayesiana e il controUo adattativo e distribuito". 
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We are interested in two problems that may arise in a network running a consensus algorithm 
when only a subset of nodes is controlled by an external input or measured by an external 
processor. Namely, is it possible to reach all the node configurations just controlling a limited 
number of nodes? Respectively, is it possible to reconstruct the entire network state just knowing 
the state of a limited number of nodes? 

In this paper we will concentrate on a first-order network system running a Laplacian based 
average consensus algorithm with fixed communication graph topology of path or cycle type. 
Average consensus has been widely studied in the last years. Several distributed feedback laws 
have been proposed. A survey on these algorithms and their performance may be found e.g. 
in [|2l and references therein. The dynamical system arising from a consensus network with 
fixed topology is a linear time-invariant system. The problem of understanding if all the node 
configurations can be reached by controlling a given subset of nodes is a reachability problem. 
Respectively, the problem of understanding if the entire network state may be reconstructed is 
an observability problem. 

We organize the relevant literature in two areas. First, the reachability problem in a first-order 
network arises in network systems where all (or most of) the nodes run a linear average consensus 
algorithm, while a subset of them can be driven by an exogenous input. These networks are 
often called leader-follower networks in the sense that the nodes controlled by an external input 
are leaders that drive the followers to desired configurations. The reachability (controllability) 
problem for a leader-follower network was introduced in ^ for a single control node. Intensive 
simulations were provided showing that it is "unlikely" for a Laplacian based consensus network 
to be completely controllable. In [4J and [l5l|, necessary conditions for controllability, based on 
suitable properties of the graph, have been provided. The conditions rely on algebraic graph 
tools based on the notion of equitable partitions of a graph. A more exhaustive analysis of 
the controllability and other structural system properties for network systems on the basis of 
graph structural properties can be found in In [|71 preliminary results are given for the 
controllability of first-order network with switching communication topology applied to the 
formation control problem. In JSl sufficient conditions for controllability in first and second order 
multi-agent systems with delay are given. The conditions are based on the eigenstructure of the 
network system and delay matrices. In ||9l necessary conditions and sufficient conditions for the 
controllability of tree graphs are given, based on the eigenvalues of suitable subgraphs. A first 
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contribution to the controllability of multi-agent systems with nonlinear dynamics is provided in 
ifTOll . where the controllability of pairs of identical nonholonomic vehicles maintaining a constant 
distance is studied. 

Second, observability for a network system running an average consensus algorithm has been 
studied for the first time in [fTTI . In that paper the authors provide necessary conditions for 
observability. The conditions are based on equitable partitions of a graph as in the reachability 
setting investigated in and [|5l. A recent reference on observability for network dynamic 
systems is [|T2l . Here, the linear dynamical systems of the network are decoupled and the coupling 
among the systems appears through the output. A parallel research line investigates a slightly 
different property called structural observability [13]. Here, the objective is to choose the nonzero 
entries of the consensus matrix (i.e. the state matrix of the resulting network system) in order 
to obtain observability from a given set of nodes. However, in many contexts the structure of 
the system matrix is given (e.g. the Laplacian for average consensus). 

It is worth noting that, the observability property is an important property in distributed 
estimation, flTH, [jTSl. and intrusion detection problems [fT6l . [fTTI for steady state analysis. 
In the literature this property is often assumed or considered as non generic. Finally, preliminary 
results on the controllability and observability of path and cycle graphs were given in fW\, 
where, using these properties, a formation control strategy was proposed. 

The contribution of the paper is twofold. First, we provide necessary and sufficient condi- 
tions based on simple algebraic relations from number theory that completely characterize the 
reachability (observability) of path and cycle graphs. More in detail, on the basis of the node 
labels and the total number of nodes in the graph we are able to: (i) identify all and only the 
reachable (observable) nodes of the graph, (ii) say if the graph is reachable (observable) from a 
given set of nodes and (iii) construct a set of control (observation) nodes from which the graph 
is reachable (observable). Interesting immediate corollaries of our results are: (i) a path graph 
is reachable (observable) from any single node if and only if the number of nodes of the graph 
is a power of two, n = 2' , i G N, (ii) a path graph is reachable (observable) from a single node 
i if and only if there is no odd prime factor p oi n such that [n — i) = [i — \) -[- ap for some 
integer a; (iii) a cycle graph is reachable (observable) from any pair of nodes if and only if n 
is a prime number, and (iv) a cycle graph is reachable (observable) from two nodes, say ii and 
Z2, if and only if Z2 — H and n + ii — i2 are coprime. Thus, e.g., any cycle is observable from 
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two adjacent nodes. 

Second, we provide a closed form expression for the unreachable (unobservable) eigenvalues 
and eigenvectors, and characterize the orthogonal complement to the reachable subspace (respec- 
tively the unobservable subspace) for any unreachable (unobservable) set of nodes. This result 
is based on the complete characterization of the spectrum of suitable submatrices of the path 
and cycle Laplacians. As a consequence of these linear algebra results, we also provide a closed 
form for all the Laplacian eigevalues of a path graph. At the best of our knowledge both the 
characterization of the Laplacian eigenvalues and the mathematical tools used to characterize 
them are new. 

The paper is organized as follows. In Section |Il] we introduce some preliminary definitions 
and properties of undirected graphs, describe the network model used in the paper and set up the 



reachability and observability problems. In Section III we provide a complete characterization 
of the eigenvalues and eigenvectors of the Laplacian of a path graph and other matrices used 
to study the path and cycle reachability (observability). Finally, in Section |IV] we provide a 
complete characterization of the path and cycle graph reachability (observability), and provide 
some useful example explaining the main results. 

Notation: We let N, Nq, and ]R>o denote the natural numbers, the non-negative integer 
numbers, and the non-negative real numbers, respectively. We denote 0^, G N, the vector 
of dimension d with zero components and 0rf^xd2' '^i;'^2 ^ the matrix with di rows and 
d2 columns with zero entries. For i G N we let be the i-ih. element of the canonical basis, 
e.g. ei = [1 . . . 0]-^. For a matrix A G M'^i^'^^ denote the {i,j)th element and 
[A]i the ith column of A. For a vector v E R.'^ we denote {v)i the ith component of v so 
that V = [(f )i . . . (f)d]^. Also, we denote 11 G R.'^^'^ the permutation matrix reversing all the 
components of v so that Uv = [{v)d ■ ■ ■ (w)i]^ (the j-th column of 11 is [n]j = e„_j+i). Adopting 
the usual terminology of number theory, we will say that k divides a nonzero integer m (written 
k\m) if there is an integer q with the property that m = kq. When this relation holds, k is said 
a factor or divisor of m. If two integers b and c satisfy for a given m the relation m\{b — c) 
then we say that b is congruent to c modulo m (written b = c mod(m) or equivalently 6™=™c). 
The greatest common divisor of two positive integers a and b is the largest divisor common to 
a and b, and we will denote it GCD{a,b). The greatest common divisor can also be defined 
for three or more positive integers as the largest divisor shared by all of them. Two or more 
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positive integers that have greatest common divisor 1 are said relatively prime or coprime. A 
prime number is a positive integer that has no positive integer divisors other than 1 and itself. 
Every natural number n admits a prime factorization (Fundamental Theorem of Arithmetic), i.e. 
we can factorize n as n = 2"° nL=i Po^' where no E N and each is an odd prime number. 
Notice that in our factorization we allow two or more factors Pa to be equal. 



In this section we present some preliminary terminology on graph theory, introduce the network 
model, set up the reachabiUty and observability problems, and provide some standard results on 
reachability (observability) of linear systems that will be useful to prove the main results of the 
paper. 

A. Preliminaries on graph theory 

Let G = (/, E) be a static undirected graph with set of nodes I = {1, . . . , n} and set of edges 
E G I X I. We denote J\fi the set of neighbors of agent i, that is, J\fi — {j e I \ e E}, 
and di — 'Ylijo^Ti ^ degree of node i. The maximum degree of the graph is defined as 
A = maxjg/ dj. The degree matrix D of the graph G is the diagonal matrix defined as [D]^ — dj. 
The adjacency matrix A e R"^" associated to the graph G is defined as 



The Laplacian L of G is defined a^ L = D—A. The Laplacian is a symmetric positive semidefinite 
matrix with k eigenvalues in 0, where k is the number of connected components of G. If the 
graph is connected the eigenvector associated to the eigenvalue is the vector 1 = [1 ... 1]^. 

We introduce the two undirected graphs that will be of interest in the rest of the paper, 
namely the path and cycle graphs. A path graph is a graph in which there are only nodes of 
degree two except for two nodes of degree one. The nodes of degree one are called external 
nodes, while the other are called internal nodes. From now on, without loss of generality, we 
will label the external nodes with 1 and n, and the internal nodes so that the edge set is 
E = + \) I i G {1, . . . , n — 1}}. A cycle graph is a graph in which all the nodes have 
degree two. From now on, without loss of generality, we will label the nodes so that the edge 
set is = {(i,imod(n) + 1) | i e {1, . . . ,n}}. 



II. Preliminaries and problem set-up 
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B. Network of agents running average consensus 

We consider a collection of agents labeled by a set of identifiers / = {!,. ..,«}, where n E N 
is the number of agents. We assume that the agents communicate according to a time-invariant 
undirected communication graph G = (/, E), where E = E I x I \ i and j communicate}. 

That is, we assume that the communication between any two agents is bi-directional. The agents 
run a consensus algorithm based on a Laplacian control law (see e.g. [2J for a survey). The 
dynamics of the agents evolve in continuous time (t E M>o) and are given by 

a^i(^) = - (a^iW - i e {1, . . . , n}. 

Using a compact notation the dynamics may be written as 

x{t) = -Lx{t), t E M>o, 

where x = [{x)i . . . = [xi . . . Xn]'^ is the vector of the agents' states and L is the graph 

Laplacian. 

Remark 2.1 (Descrete time system): In discrete time, we can consider the following dynamics 

Xi(t + 1) = - e ^ (x,;(t) - Xj(t)), i E {1, . . . ,n}, 

where e G M is a given parameter. A compact expression for the dynamics is 

x{t +!) = (/- eL)x{t), t E No. 

For e G (0, 1/A) (A is the maximum degree of the graph), P=(I — eL) is a nonnegative, doubly 
stochastic, stable matrix. 

It can be easily shown that the continuous and discrete time systems have the same reachability 
and observability properties (namely the same unreachable and unobservable eigenvalues and 
eigenvectors). Therefore, the results shown in the paper also hold in this discrete time set-up. □ 



C. Network reachability and observability 

In this section we describe the mathematical framework that we will use to study the reacha- 
bility and observability of a network system. We start by describing the scenarios that motivate 
our work. As regards the reachability problem, we imagine that in a network of agents running 



average consensus as in Section II-B[ a subset of nodes can be controlled by an external input. 
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In the literature these nodes are often called leader nodes or pinned nodes. The idea is that they 
are special nodes with higher computation capabilities so that more sophisticated control laws 
can be designed. Let = {ii, . . . , C {1, . . . , n} be the set of control nodes, the dynamical 
system modeling this scenario is 

x{t) = -Lx{t) + Bu{t), 

where u{t) = . . .,Ui^{t)f is the input and B = [d^ | . . . | dj. 

As regards the observability problem, we imagine that an external processor (not running the 
consensus algorithm) collects information from some nodes in the network. We call these nodes 
observation nodes. In particular, we assume that the external processor may read the state of 
each observation node. Equivalently, we can think of one or more observation nodes, running 
the consensus algorithm, that have to reconstruct the state of the network by processing only 
their own state. We can model these two scenarios with the following mathematical framework. 
For each observation node i G /, we have the following output 

yi{t) = x,{t). 

Given the set of observation nodes Jq = {ii,...,im} C {l,...,n}, the output is y(t) = 

1 ^ 

i't) i't) ■ ■ ■ (t) . Therefore, the system dynamics is given by 

x(t) = -Lx{t), 

y{t) = Cx{t), 
where the output matrix is C = [cj^ | • • • | Cj^]^. 

It is a well known result in linear systems theory that the reachability properties of the pair 
{L,B) are the same as the observability properties of the pair {L,C) = = (L,B'^). 

Remark 2.2 (Duality and regulator design): Due to the symmetry of the state matrix L, the 
network is reachable from a given subset of nodes if and only if it is observable from it. This 
important property allows, for example, to design a regulator at a subset of (control/observation) 
nodes that estimates the entire network state and controls it to a desired configuration. This has, 
for example, certainly an impact on security issues. □ 

Remark 2.3 (Equivalence with other problem set-ups): Straightforward results from linear sys- 
tem theory can be used to prove that the controllability problem studied, e.g., in [4] and [5] and 
the dual observability problem studied in [[TT]| can be equivalently formulated in our set up. □ 



September 19, 2011 



DRAFT 



8 



D. Standard results on reachability and observability of linear systems 

The reachability problem consists of looking for those states that can be reached in finite time 
from the origin. Respectively, the observability problem consists of looking for nonzero values 
of x{Q) that produce an identically zero output y{t). Using known results in linear system theory 
the two problems are equivalent to studying the rank of the reachability matrix, 

nn=[B\LB\ ... 

respectively of the observability matrix 

C 
CL 



- n—l 



On 



CL 



n-l 



The image Xr of TZn is the reachable subspace, i.e. the set of states that are reachable from 
the origin. The kernel X„o of On is the unobservable subspace, i.e. the set of initial states that 
produce an identically zero output. 

Here, we recall an interesting result on the reachability (observability) of time-invariant linear 
systems known as Popov-Belevich-Hautus (PBH) lemma. 

Lemma 2.4 (PBH lemma): Let A e R"""", B e W""" and C e R"*^", n,meN, be the state, 
input and output matrices of a linear time-invariant system. The pair {A,B) is reachable if and 
only if 



rank 



A-XI\B 



for all A e C. Respectively the pair {A, C) is observable, if and only if 

C 



rank 



A -XI 

for all A e C. □ 
Combining the PBH lemma with the fact that the state matrix is symmetric (and therefore 

diagonalizable) the following corollary may be proven. 

Corollary 2.5 (PBH lemma for symmetric matrices): Let A G W""'', B e R"^™ and C e 

j^mxn^ n, m e N, be the state, input and output matrices of a linear time-invariant system, where 
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A is symmetric. Then, the orthogonal complement to the reachable subspace X^, associated to 
the pair {A,B), is spanned by vectors vi satisfying 

(1) 

Lvi = \vi. 

Respectively, the unobservable subspace X„„ associated to the pair (L, C) is spanned by vectors 
vi satisfying 

CVI = Om 

(2) 

Lvi = Xvi, 

for A G M. □ 
In the rest of the paper we will denote the eigenvalues and eigenvectors for which ([!]) is 
satisfied, unreachable eigenvalues and eigenvectors. Respectively, we will denote unobservable 
eigenvalues and eigenvectors the ones for which (|2]) is satisfied. 

III. Spectral properties of the Laplacian of a path and related submatrices 

In this section we provide a closed form expression for the eigenvalues and eigenvectors of 
suitable submatrices of the Laplacian of path and cycle graphs. This characterization will play 
a key role in the characterization of the reachability (observability) properties of path and cycle 
graphs. As a self-contained result of this section, we provide a closed form expression for the 
Laplacian spectrum of a path graph. 

We start motivating the analysis in this section. Recall that, without loss of generality, we 
assume that nodes of the path are labeled so that the undirected edges are {i, i + 1) for i G 
{l,...,n — 1}. Let Ln denote the Laplacian of a path graph of length n and B (C) the input 
(output) matrix associated to the set of control (observation) nodes Jq = {ii, . . . , v}- Using the 



PBH Lemma in the version of Corollary 2.5[ unreachability (unobservability) for the path graph 



from lo is equivalent to the existence of a nonzero solution of the linear (algebraic) system 
LnV = \v, where v satisfies B'^ v = {C v = 0), i.e. {v)j = for each j G /„. Exploiting the 
structure of the linear system we can write 
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1 -1 ... 

-1 2 

-1 

OH 









-1 
1 



" (^')l " 






(^')n-i 




(^)n-i 





= A 





(^^)ii+i 






. (^)" . 







t ii-th column 

where the vertical line on the zi-th column means that it is multiplied by {v)i-^ = 0. The same 
holds for each j-th column, j e Ig. Now, define the matrices N,, e W^'^ and e M'^^^ as 



1 -1 

-1 2 

■•• -1 

-12 



and M,, 



2-1 

-1 2 

■•• -1 

-12 



where the subindex refers to their dimensions. The Laplacian L„ can be compactly written as 

J) ... 













-^i2-w-l 









nNn_i n 



■il-th column ^ (12 — ii)-th column 'j^ . . . 

where 11 = 11^ = 11" ^ is the (symmetric) permutation matrix reversing all the components of a 
vector. 
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Remark 3.1 (Partition of the Laplacian of a cycle): Applying the same procedure to the Lapla- 
cian of a cycle, under the agreement of labeling the nodes so that ii = 1, we get a partition of 
the Laplacian where the submatrices are all matrices of type M^, fi eN. □ 

We are now ready to investigate the spectral properties of these matrices. We begin with a 
useful lemma. 

Lemma 3.2: The eigenvectors of A*"^, and L^, /i G N, have nonzero first and last compo- 
nents. 

Proof: We prove the statement by proving that if an eigenvector has zero first component, 
then all the other components have to be zero. The same line of proof can be followed to prove 
that if the last component is zero, then the same must hold for the previous /i — 1. 

Let A be an eigenvalue of A^"^ (respectively and L^J and v = [vi . . . v^] the corresponding 
eigenvector. The following conditions must hold 

avi — V2 = Af 1 

-Vi_i + hvi - Vi+i = Xvi, « G {2, . . . , /i - 1} 

-v^,-^i + bv^ = XVf,, 

for suitable (positive) a and b. We proceed by induction. We use the inductive assumption that 
Vj^i = Vj = and prove that vj^i = 0. The statement is obviously true for j = 1. Indeed, from 
the first equation it follows easily that if i;i = then f 2 = (we have considered a fake vq = 
in the first equation). Then, plugging the inductive assumption in the second equation, we get 
Vj^i = for j G {2, . . . , /i — 1}, while the last equation gives the result for j = /i. ■ 
Remark 3.3 (A path is reachable (observable) from an external node): Combining the previ- 



ous lemma with Corollary 2.5 , it follows easily that a path graph is reachable (observable) from 
each of the external nodes as shown, e.g., in flU ( [[T9l ). □ 
Proposition 3.4 (Eigenstructure of and M^): For any two matrices N^, G W^" and G 
^MXM the following holds: 

(i) All the eigenvalues of Ny are eigenvalues of and the corresponding eigenvectors, 
respectively f G M'^ and w G M , are related by w 
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(ii) Eigenvalues and eigenvectors of Ny and have the following closed form expression: 



2-2COS [(2/c-l)^] 
(z/ + j)(2A;-l)7r 



sm 



2z/ + 1 



(3) 



2-2COS (fc^) 
jkiT 



sm 



/i + 1 



A; 



(4) 



Wi 
W2 



V 



with V G 



Proof: To prove the first statement, consider a vector w = 

We show that u' is an eigenvector of with eigenvalue A at if and only if f is an eigenvector 
of Ni, with eigenvalue X^. Indeed, by exploiting the structure of we get 

Ny + eief 



2u 



Wi 
W2 



T 



Uv 

V 



UNyU Uv + e^ie^U v - ej v) 
-ei(e^n v) + NyV + ei{v)i 



UNyV + ey[{v), - iv)i] 




UNyV 


Nyv + ei[{v)i - {v)i] 




Nyv 



Now, notice that, if A at is an eigenvalue of N^, i.e. N^v = Xnv for some nonzero v G W, then 

Uv 



w is an eigenvector of M2u with the same eigenvalue. Viceversa, if w 



7^ satisfies 
lus concluding 



M2uW = Xnw, then, from the last u relations, it follows easily that N^^v = X^v, t 
the first part of the proof. 

To prove the second statement, we observe that = 2/ — Ap, where Ap is the adjacency 
matrix of the path, reported in Appendix |Aj Thus, the eigenvalues of A<f^ can be computed by 
summing 2 to the eigenvalues of Ap and the eigenvectors are the same, so that Q follows. We 
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now show that the v eigenvalues and eigenvectors of are the eigenvalues and eigenvectors 
of corresponding to the odd values of k in (|4]). To see this, notice that the v eigenvectors 

of with structure w = satisfy {w)j = {w)2v-j+i, j = l,...,z/. Referring to the 

V 

parametrization of w in Q, we look for ks satisfying sin (^^^) = sin (^ ''^'^2i^^i''^"^ )' ^ 
{1, . . . ,1^}. This is easily verified for odd ks because the arguments of the sine functions on the 
two sides sum to (21 + l)n for some integer /. This concludes the proof. ■ 

Next, the eigenvalues of the Laplacian L„ of a path graph of length n are expressed in closed 
form by relating them to the eigenvalues of the M„ i matrix. The following technical lemma 
gives the tools to directly compute the eigenvalues of the laplacian matrix of a path graph. 

Lemma 3.5: The characteristic polynomials of A^"^ and M^, yU G N, satisfy the following 
relations, 



det(s/ - 


N,) = 


(s - 1) det(s/ - 


M^_0 


- det(s/ - 


M^_2 


det(s/ - 


N,) = 


[s - 2) det(s/ - 


N,-i) 


- det(sJ - 




det(s/ - 


M,) = 


[s - 2) det(s/ - 




- det(s/ - 


M^_2 



Proof: The result is obtained by applying the Laplace expansion to compute the determinant 
of the matrices (si — N^j) and (si — M^). In particular, the first and second relations are obtained 
by expanding respectively along the first and last row of (si — N^j). The third relation is obtained 
by expanding along, e.g., the first row of (si — M^). ■ 
Proposition 3.6 (Eigenvalues of the Laplacian of a path): The characteristic polynomial of the 
Laplacian, L„, of a path graph of length n can be written as 

det(sJ - L„) = s det(s/ - M„_i). 

Thus, the eigenvalues of the Laplacian are given by 

Al„ = 2-2cos(^(A;-l)^) ,fc = 

Proof: Applying the Laplace expansion for the computation of the determinant and using 
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the results of the previous lemma, the following equalities hold: 

det(s/ - L„) = (s - 1) det(s/ - iV„_i) - det(s/ - Nn-2) 

= {s- if det(s/ - Mn-2) - 2(s - 1) det(s/ - M„_3) - det(s/ - M„_4) 
= (s^ - 2s) det(s/ - Mn-2) - sdet(s/ - M„_3) 

= (s^ - 2s) det(s/ - Mn-2) - s [(s - 2) det(s/ - Mn-2) - det(s/ - M„_i)] 
= sdet(s/-M„_i). 

■ 

IV. Reachability and observability of path and cycle graphs 

In this section we completely characterize the reachability and observability of path and cycle 
graphs. 

A. Reachability and observability of path graphs 

We characterize the reachability (observability) of a path graph by using the PBH lemma in 
the form expressed in Corollary 2.5[ First, as recalled in Remark 3.3[ a path graph is always 



(reachable) observable from nodes 1 or n. Next two lemmas give necessary and sufficient 
conditions for reachability (observability) from a given subset of nodes in terms of the Ny 
and M^ submatrices introduced in the previous section. 

Lemma 4.1: A path graph of length n is reachable (observable) from a node i E {2, . . . ,n—l} 
if and only if the matrices Ni-i and Nn-i do not have any common eigenvalue. The eigenvalues 
common to the two matrices are all and only the unreachable (unobservable) eigenvalues of the 
Laplacian L„ from node i. 



Proof: Applying Corollary 2.5, we have that i is not reachable (observable) if and only 
if LnV = \v and efv = for some A and v. Equivalently, i is not reachable (observable) if 
and only if there exists an eigenvector v of L„ with v = [vi f2]^, fi G W^^, V2 G M"^^. 
Component-wise this is written as 

Ni_iVi = Xvi 

{vi)i-i + {V2)i = (5) 

{UNn-iU) V2 = XV2. 
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The necessary condition follows easily by the above equations. Indeed, if Ni_i and Nn-i have 
at least one common eigenvalue Aq with corresponding eigenvectors respectively fio and V20, 
then the conditions in the above equations are satisfied for v = [vio pv2o] and A = Aq, where 
p G M is just a scaling factor to satisfy (fio)i-i + p('i^2o)i = 0. 

To prove the converse we proceed by analyzing when these three conditions are satisfied. 
The first equation in Q is verified in two cases: i) f 1 = and A arbitrary, and ii) vi and 
A respectively eigenvector and eigenvalue of A^i-i. From the first condition it follows easily 



that (^2)1 = and, using Lemma 3.2. that t>2 = 0. Therefore, the only possibility to have 
unreachability (unobservability) is ii). Now, using the second equation it follows easily that 
V2 7^ 0, and, using the third equation, that Ilv2 must be an eigenvector of Nn-i corresponding 
to the same eigenvalue A of Ni_i. This concludes the first part of the proof. 

The fact that the unreachable (unobservable) eigenvalues of L„ from node i are all and only 
the eigenvalues common to A^j_i and Nn-i follows straight from the previous argument. Indeed, 
by definition the unreachable (unobservable) eigenvalues are all and only the ones satisfying the 
condition in Corollary 2.5 and, thus, the equations in ([5]). ■ 

Remark 4.2 (Paths with odd number of nodes): A straightforward consequence of the previ- 
ous lemma is that a path graph with an odd number, n, of nodes is not reachable (observable) 
from the central node. Also, {n — l)/2 eigenvalues are unreachable (unobservable) from that 
node, namely the eigenvalues of L„ that are also eigenvalues of A''(„-i)/2- The corresponding 
{n — l)/2 unreachable (unobservable) eigenvectors of L„ are of the form [v^ -t>^n]^ with v 
being the eigenvectors of Ni^n~i)/2 as in ([3]). □ 

A generalization to the multi input (output) case is given in the following lemma. 

Lemma 4.3: A path graph of length n is reachable (observable) from the set of control (ob- 
servation) nodes lo = {ii, . . . ,im.} if and only if the matrices iVj^_i, Mi^-i^^i, . . ., Mi„__i^_^_i 
and Nn-i^ do not have common eigenvalues. The eigenvalues common to the matrices are all 
and only the (unreachable) unobservable eigenvalues of L„ from the set I^. 

Proof: We proceed as in the single input (output) case and apply the PBH Lemma in 



the version expressed in Corollary 2.5 The path is not reachable (observable) from the set 
lo = {^1, • • • 5 irn] if and only if there exists an eigenvector v of L„ with {v)i^ = . . . = (f = 0, 
so that f = [ff fj... f^]^ for suitable vectors vi G M*^"-*^, . . . , e IR*™"^ This is 
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equivalent to 



Ni^^iVi = Xvi 







K)n-l + (^^2)1 

{V2)i2-1 + (^3)1 = 



(6) 



(^^m-l)j„-l + {Vm)l = 



The proof follows by using the same arguments in Lemma 4.1 



We are now ready to completely characterize the reachability (observability) of a path by means 
of simple algebraic rules from number theory. For the sake of clarity, we state the theorem for 
path graphs of length n, where n has a prime factorization with distinct odd prime factors. The 
general case follows straight and is discussed in a remark. 

Theorem 4.4 (Path reachability and observability): Given a path graph of length n, let n = 
2"o 11^=1 Pi' be a prime number factorization for some /c G N and distinct (odd) prime numbers 
Pi, . . . ,Pk. The following statements hold: 

(i) the path is not completely reachable (observable) from a node i G {2, ... ,n — 1} if and 
only if 

/ -x modp , . . 

[n — l) = [t — 1) 

for some odd prime p dividing n; 

(ii) the path is not completely reachable (observable) from a set of nodes lo = {n, 
{2, . . . , n — 1} if and only if 



2{h - 1) + 1 



mod p , 



1^2 - 



, modp modp. 



modp . ^ . 

^m— 1 zyn 



+ 1, 



for some odd prime p dividing n; 
(iii) for each odd prime factor p G {pi, . . . ,pk} of n, the path is not reachable (observable) from 
each set of nodes = {£p— ^2^}£e{i,...,-} with the following unreachable (unobservable) 
eigenvalues 



= 2 - 2 cos (2z/ - 1 



p 



z/G{l,...,^h— }; 



(7) 
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and unreachable (unobservable) eigenvectors 



(8) 



vl -{Uv.f -vl ... {-ly^iliv,) 

where the eigenvector of A''(p-i)/2 corresponding to the eigenvalue \y for 

z/G{l,...,(p-l)/2};and 
(iv) if node i belongs to ll^ = {igj — ... jl} for / < k distinct prime factors qi ^ 

... ^ qi of n, then the set of unreachable (unobservable) eigenvalues from node i is given 
by 

A. = 2-2cosf(2z.-l) ^ V z.G{l,...,^^^-^^-*l^}. 



qi- ... -qij 2 
Also, the orthogonal complement to the reachable subspace, (X^)-*-, (respectively the 
unobservable subspace, Xno) is spanned by aU the corresponding eigenvectors of the form 

1 T 



7j -(n^;,)^ ... {-l)HUv,r 



where y,^ G ]R((''i is the eigenvector of A^((gj.....g;)_i)/2 corresponding to the eigen- 
value Aj, for G {1, . . . , ((gi ■ . . . ■ g,) - l)/2}. 
Proof: Using Lemma |4]T we have that the path graph is not completely reachable (observ- 



able) from node i if and only if Ni_i and N„_i have at least one common eigenvalue. Therefore, 
using Proposition 



3.4 



we have that it must hold 2 — 2 cos = 2 — 2 cos |?^^-4tt, for 

2(2—1)4-1 2(n— jj+l' 

some ji G {1, . . . ,i — 1} and j2 ^ {1, • • • , — i)}- In the admissible range of ji and j2 the 
cosine arguments are less than vr so that the cosine is one to one. Thus, the equality holds if 
and only if 2^_^i)+i ~ ■y^n-^yli ■ integers ji and j2 satisfying this equation exist in the 
admissible range if only if 2(i — 1) + 1 and 2(n — i) + 1 are not coprime, that is, if and only if 
GCD{2i - 1, 2n - 2i + 1) is greater than one. Now, GCD{2i -l,2n-2i + 1) is odd because 
2i — 1 and 2n — 2i + 1 are. Therefore, we can write 2i — 1 = pai and 2n — 2i + 1 = pa2 with 
p, Ci\ and odd. This is equivalent ton + ra — 2z + l= poL2 and, since p divides ra, p must 
divide also (n — 2i + 1) thus concluding the first part of the proof. 



To prove statement (ii), we have by Lemma 4.3 that the path graph is not reachable (observable) 



from the set = {ii, . . . , z„} if and only if the matrices Ni^_x, Mja-n-i' • • and 



Nn~i^ do not have common eigenvalues. The proof follows by using again Proposition 3.4 and 
the same arguments as in the single node case. 

September 19, 2011 DRAFT 



18 



To prove statement (iii), we start observing that the set of nodes = {£p — ^2^}^g{i ... n}, is 



the set of all nodes satisfying condition in (i) for a given p G {pi, . . . ,Pk}- Using Lemma 4.3 



we have that the unreachable (unobservable) eigenvalues from this set of nodes are the common 



eigenvalues to iV(p_i)/2 and Mp_i. From Proposition 3.4 it follows easily that the common 
eigenvalues between A''(p_i)/2 and Mp are all the eigenvalues of A^(p_i)/2 and have the form 
in equation As regards the unreachable (unobservable) eigenvectors, using equation ([6]) 
for this special set I^, it follows straights that the eigenvectors have zero components as in 
equation ([8]). To prove that the nonzero components have that special structure in equation ([8]), we 
observe that they must be eigenvectors of respectively N(^p^iy2, Mp_i, . . .,Mp_i and niV(p_i)/2n. 



From point (i) of Proposition 3.4, the eigenvectors of A''(p-i)/2, Mp_i, . . .,Mp_i and nA''(p_i)/2n 
are respectively aoVi^, ai[{Ilv,y)^ "^lY ^ • • OLik_i[{JlVyf v^Y ^i^d anliv^, where Vy is an 
eigenvector of N(^p^i)/2 and «^ G M, /U G {0, . . . , ^}. Finally, using again equation ([6]), = 
— q;^_i, /iG{l,...,^}, so that the proof follows by choosing = 1. 

The proof of statement (iv) follows from the definition of unreachable (unobservable) eigen- 
values and eigenvectors, and arguments as in the previous statements. ■ 



4.4 



Remark 4.5 (General version of Theorem 4.4): In the general case of a path graph of length 
n = 2^° Y[t=iP'^' where pi, . . . ,pk are not all distinct, statement (i) and (ii) of Theorem 
continue to hold in the same form. As regards statement (iii), it still holds in the same form, 
but it can also be strengthen with a slight modification. That is, for each multiple factor p 
with multiplicity k, the statement continues to hold if p is replaced by with a G {1, . . . ,k}. 
Statement (iv) holds if for each prime factor p with multiplicity k we check if node i belongs not 
only to but also to each I^" with a G {1, . . . , A;}. Consistently the unreachable (unobservable) 
eigenvalues and eigenvectors considered in the statement must be constructed by using p" instead 
of p, where d = maxa{a G {1, . . . , A;}|i G /^"}. □ 



The following corollary follows straight from Theorem 4.4 and characterizes all and only the 
path graphs that are observable from any node. 

Corollary 4.6 (Reachable (observable) paths from any node): Given a path graph of length 
n = 2^ for some A; G N, then the path is reachable (observable) from any node. □ 

Next, we provide a simple routine giving a graphical interpretation of the results of the theorem. 
We describe the routine for paths with simple factorization leaving the generalization to a specific 
example. We proceed by associating a unique symbol to each set of nodes defined at point 
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(iii) of Theorem 4.4 for a given p. Thus, each group of nodes sharing the same unreachable 
(unobservable) eigenvalues has the same symbol and nodes of different groups (and so associated 
to different unreachable (unobservable) eigenvalues) have different symbols. Formally, let n = 
2"o Y[t=iPf some no G N and pi, . . . ,pk prime integers. At the beginning of the procedure 
we initialize all the nodes without any symbol. For any p^, v G {1, . . . , A;}, we partition the 
nodes into njpy groups of py nodes and assign the same symbol to all the nodes in position 

i = jpi, — ^^Y^, j G {1, . . . , ^}. A set of nodes from which the path is reachable (observable) 
is obtained by selecting any node without symbols, if there are any, or a set of nodes having no 
symbols in common. 

Three examples for n = 6 (even), n = 15 (odd) and n = 9 (multiple factor) are shown 
respectively in Figure [Tj in Figure [2] and in Figure |3] In Figure [T] nodes with the triangle symbol 
are unable to reconstruct the state of the network by themselves. Indeed, they share the same 
unreachable (unobservable) eigenvalue A = 1. In view of the previous results, focusing on node 

ii = 2, notice that A^j_i = Ni = [1] (whose eigenvalue is 1), A^„_i = A^4 and its eigenvalues 
are: {0.12,1,2.35,3.53}. The common eigenvalue is of course A = 1. The unreachability 
(unobservability) can be checked more easily using the test {n — i) = 4 ™= 1 = {i — 1). 



Fig. 1. Reachable (observable) nodes for a path with n = 6 nodes. 



In Figure [2] nodes with the triangle belong to the set J^^, with pi = 3, and nodes with the 
square to I^^, with P2 = 5. The "triangle nodes" share the same (unreachable) unobservable 
eigenvalue A = 1 (A^j-i = Ni = [1]), while two unreachable (unobservable) eigenvalues, 0.3820 
and 2.6180, are associated to the "square nodes". Finally, the central node has both the triangle 
and square symbols. 

In Figure |3] we consider a path of length n = 9 = 3^, with p = 3 being a multiple factor. 
Nodes with the triangle belong to the set I^^, with pi = 3. To the central node is associated both 

2 

a triangle (since it belongs to /^) and a square since it is the unique node in /o\ with pf = 9. The 
"triangle nodes" share the same (unreachable) unobservable eigenvalue A = 1 (A^j_i = Ni = [1]), 
while the central node has four unreachable (unobservable) eigenvalues being the eigenvalues of 
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Fig. 2. Reachable (observable) nodes for a path with n = 15 nodes. 

A^4. This example suggests how to associate symbols in the general case when n has multiple 
factors. That is, for each multiple factor p with multiplicity k, we use k different symbols, one 
for each p"', a E {I, . . . , k}. 

1 23456789 

A AD A 

Fig. 3. Reachable (observable) nodes for a path with n = 9 nodes. 

B. Reachability and observability of cycle graphs 

Next, we characterize the reachability (observability) of a cycle graph. We start with a negative 
result, namely that a cycle graph is not reachable (observable) from a single node. First, we need 
a well known result in linear systems theory [|20| . 

Lemma 4.7: If a state matrix A G M"^", n G N, has an eigenvalue with geometric multiplicity 

> 2, then for any i? G M" (respectively C G M^^") the pair (A, B) is unreachable (respectively 
the pair {A, C) is unobservable). □ 

As shown in Appendix |B} all but at most two eigenvalues of the Laplacian of the cycle have 
geometric multiplicity two. Thus, applying the previous lemma next proposition follows. 

Proposition 4.8: A cycle graph is not completely reachable (observable) from a single node 
for any choice of the control (observation) node. Furthermore, \^^~\ eigenvalues are unreachable 
(unobservable). 
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Proof: The unreachability (unobservability) follows directly by noting that all but at most 
two eigenvalues of the cycle have geometric multiplicity two and by applying the previous lemma. 
By using the PBH lemma, it follows straight that the unreachable (unobservable) eigenvalues 
are all but the zero eigenvalue. ■ 

Recall that, without loss of generality, we label the nodes of the cycle so that the undirected 
edges are [i mod (n) + 1)) for i E {1, . . . ,n}. Following the same line as in the reachability 
(observability) analysis of path graphs, we can prove the following lemma. 

Lemma 4.9: A cycle graph of length n is reachable (observable) from the set of nodes lo = 
{ii, ...,im} if and only if the matrices Mi^^i^^i, . . ., Mj^_i^_^_i and M(i^_i^_i)modn do not 
have common eigenvalues. 

Proof: Without loss of generality, we set node ii = 1, so that the M matrices in the 
statement of the theorem become Mi^-i-^-i, . . ., Mj^_j^_^_i and Mn-i^. Following the same 
lines of the proof of Lemma |43 loss of reachability (observability) is equivalent to the existence 



of a nonzero solution f = [0 f J . . . f^]^ to: 

(fl)l + (t^m)n-i„ = 

Mi2_i^_lfl = Xvi 



Now, if Vj = for some j G {1, . . . , m}, using Lemma 3.2 the only solution of the above system 
is f = 0. With this condition in hand, it follows straight that t> is a (nonzero) solution if and 
only if all the Vj are nonzero eigevectors of Mjj-ii-i, . . . , M^-i^ with common eigenvalue A, 
thus concluding the proof. ■ 
It is worth noting that, due to the symmetry of the cycle, the reachability (observability) 
properties are determined by the relative distance between each pair of consecutive control 



(observation) nodes. The following theorem parallels Theorem 4.4 As for the path, we state 



the theorem for cycle graphs of length n, where n has a prime factorization with distinct prime 



factors. The general case follows straight from similar arguments as in Remark 4.5 



Theorem 4.10 (Cycle reachability and observability): Given a cycle graph of length n, let 
n = Y[t=iPi' be a prime number factorization for some k E N and distinct prime numbers 
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Pi, . . . ,Pk (including the integer 2). The following statements hold: 

(i) the cycle graph is reachable (observable) from the set of nodes Jq = {ii, . . . , i^} if and 
only if 

GCD{{i2 - ii), {is - i2), ...,{n + H- ij) = 1; (9) 

(ii) for each prime factor p of n and for each fixed k E {1, . . . ,p}, the set of nodes = 
{/€ + £p}£g{o,...,ii-i}, is unreachable (unobservable) with the following unreachable (unob- 
servable) eigenvalues 



IT 

Xy = 2 — 2 cos ( V — 
p 



V e 1} 



and, for K = 1, eigenvectors 



w 



w 



T 



1 T 



(10) 



(11) 



where Wy e is the eigenvector of M(p_i) corresponding to the eigenvalue A,, 

(iii) if a set of control (observation) nodes Jq with cardinality greater than 1, satisfies 



GCD[{i2 - ii), (^3 - ^2), . . . , (n + ii - ira)) = qi ■ 



Qi 



where qi ^ ■ ■ ■ qi are / < k distinct prime factors of n, then the set of unreachable 
(unobservable) eigenvalues from lo is given by 



Xy = 2 — 2 cos ( u 



I, y e{l,...,{qi- ...-qi)-!}. 

qi- ■■■■qij 

Also, without loss of generality, setting ii = 1, the orthogonal complement to the reachable 
subspace, (X^)-*-, (respectively the unobservable subspace, X„o) is spanned by all the 
corresponding eigenvectors 

1 T 



w. 



w. 



where the eigenvector of corresponding to the eigenvalue 

A,. 

Proof: We provide just a sketch of the proof since it follows the same line as the proof of 



Theorem 4.4 Statement (i) is proven by using Lemma 4.9 and the structure of the eigenvalues of 



the M matrices given in Proposition 3.4 with the same argument as in Theorem 4.4 (i) and (ii). 
To prove statement (ii), we start observing that the set of nodes = + ... s}, is the set 
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of all nodes satisfying condition in (i) for a given p E {pi, . . . ,Pk}- Using Lemma 4.9 we have 



that the unreachable (unobservable) eigenvalues from this set of nodes are the eigenvalues of 



Mp_i. The proof follows by the same arguments as in Theorem 4.4 (iii). Finally, statement (iii) 
follows straight. ■ 

Next corollaries provide respectively an easy way to choose two control (observation) nodes 
to get reachability (observability) for any cycle length and the class of cycle graphs (lengths) 
for which reachability (observability) is guaranteed for any pair of nodes. 

Corollary 4.11: Any cycle graph is reachable (observable) from two adjacent nodes. □ 

Corollary 4.12: A cycle graph of length n is reachable (observable) from any pair of nodes 
if and only if n is prime. □ 

As for the path, we provide a simple routine giving a graphical interpretation of the results 
of the theorem. Again, we give the procedure for the case of simple factors. 

We will mark each unreachable (unobservable) node with a different symbol. Let n = Y[t=iPi' 
for some k eN and pi, . . . ,pk distinct prime integers (here we include 2 among the p^, as well). 
At the beginning of the procedure all the nodes are initialized without any symbol. For any p^, 
v E {I, . . . , k} and i E {1, • . . ,Pu}, partition the nodes into n/pi, groups of p„ nodes and assign 
the same symbol to all the nodes in position i + k-p^, j E {1,...,^ — 1}. Nodes with the same 
symbol have the same unreachable (unobservable) eigenvectors, in the sense that controlling 
(observing from) all the nodes at the same time gives the same unreachable (unobservable) 
eigenvectors. A set of nodes from which the cycle is reachable (observable) is obtained by 
selecting any subset of nodes having no symbols in common. 

In Figure |4] there are two symbols for each node. This is because n = 15 = 3 ■ 5. Symbols 
closer to the nodes have periodicity 5 and the others have periodicity 3. Notice the ease of design 
using the above procedure. For example {4, 13} and {8, 14} are unreachable (unobservable) pairs 
since they share respectively the square and the parallelogram, while {2, 13} and {5, 12} are 
reachable (observable) pairs. Finally, notice that two neighboring nodes have always different 



symbols in accordance with the result in Corollary 4.11 



V. Conclusions 

In this paper we have characterized the reachability (observability) of path and cycle graphs 
in terms of simple algebraic rules from number theory. In particular, we have shown what are 
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Fig. 4. Graphical interpretation of the observability of a cycle with 15 nodes. 



all and only the unreachable (unobservable) set of nodes and provided simple routines to select 
a set of control (observation) nodes that guarantee reachability (observabihty). 

Promising avenues for further research include the extension of the proposed methodologies 
to more complex graphs having paths and cycles as constitutive graphs (e.g., grid, torus and 
cylinder graphs). 



Appendix 

A. Adjacency matrix of path graphs 

Eigenvalues and eigenvectors of the adjacency matrix, Ap E M"^", of a path graph of length 
n can be easily computed f2T\ . Here we briefly summarize some of the steps. 
Consider the matrix 



Ar. 



1 

1 





1 







_0 ... 1 0_ 

By definition, an eigenvalue A and a corresponding eigenvector i; G M" of Ap satisfy {v)i-i — 
X{v)i + = with (f)o = (f)n+i = and at least one {v)i,i G {1, • . . ,n}, nonzero. Set 

{v)i = a and build the sequence {v)i according to (f)i+i = — (f)i-i + X{v)i (e.g. {v)2 = Aa, 
(17)3 = A^a - a, {v)4. = X^a - 2Xa, {v)i = p£{X)a = (A ■ P£-i{X) - p^_2(A))a). Imposing 
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Pn+iW = one finds = 2 cos (k^;:^) as all possible values giving (f)„+i = with nonzero 
a. Now, choose = sin (k^^), the corresponding eigenvector can be expressed componentwise 
{vk)i = sin (i ■ k^^) according the recursive formula above and simple trigonometric rules. 



B. Circulant matrices and eigenstructure of the Laplacian of a cycle graph 
An n X n matrix C of the form 

Co Cn-l ... 02 Cl 
Cl Co C„_i 02 

■ Cl Co ■ ■ : 



c 



On-2 



On-1 



On-l On-2 ■ ■ ■ Ci Cq 

is called a circulant matrix [|22| . A circulant matrix is fully specified by the first column c = 
[co, . . . , On-iY of C. The other columns are obtained by a cyclic permutation of the first. The 
eigenvalues of a circulant matrix can be expressed in terms of the coefficients Cq, . . . , c„_i [22J: 
Xj = ^^Zluj^'^Ok, ui = e*^, where here i represents the imaginary unit. 

The Laplacian matrix of a cycle graph is a special case of this family corresponding to cq = 2, 
Cl = c„_i = -1, c,- = 0, j = 2, n- 2: 



5 



A, 



yjjo^ _ e^^^>-i) - e^"^- = 2 - 2 cos f —j 

n 



3 



0, 1, ... ,n — 1. Notice that the eigenvalues Ao = and An = 4 (only if n is even) are 



simple, all the others verify Xj = X 



n-j- 
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